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I. INTRODUCTION 

The problem of a classical heavy particle embedded in 
a quantum environment is a paradigm that can be ap- 
plied to diverse physical systems. The condition for its 
applicability is the existence of a macroscopic variable 
that can be treated as classical, coupled to quantum de- 
grees of freedom. If the system allows for a separation 
of time scales such that the characteristic times of the 
quantum degrees of freedom are much faster than the 
classical ones, the evolution can be described within an 
adiabatic expansion, in which the velocity of the classical 
variable is taken as a small parameter. The Hamiltonian 
of the quantum system becomes parametrically depen- 
dent on time through the classical degrees of freedom. 
As the states of the quantum system evolve in time, they 
acquire a geometric phase, denominated Berry phase, in 
addition to the usual dynamical phased 

The backaction of the quantum environment on the 
classical degrees of freedom can be cast in terms of ef- 
fective reaction forces that affect the dynamics of the 
classical variables. The simplest and best known of these 
reaction forces is the Born-Oppenheimer force associated 
with the adiabatic potential surfaces of the fast quan- 
tum system as function of the slow classical variables. 
The Born-Oppenheimer force depends only on the coor- 
dinates of the classical degrees of freedom and is inde- 
pendent of their velocity. As emphasized by Berrjd and 
othergMS', additional reaction forces appear when going 
to next order in the adiabatic approximation, retaining 
forces which are linear in the velocity of the classical 
variables. In fact, they found that the Berry phase is 
mirrored by a Lorentz-like force, which was dubbed "ge- 
ometric magnetism." It is not associated to a real mag- 
netic field, but to an emergent geometrical property of 
the Hilbert space. Like the true magnetic Lorentz force, 
this emergent Lorentz force is non-dissipative. In general, 
one may also expect a friction force in linear order in the 




FIG. 1: Example of a scattering system that is coupled to slow 
classical degrees of freedom X(i). The movement of the scat- 
terer changes the scattering potential V[X(£)]. The backac- 
tion of the electrons passing through the scatterer then leads 
to reaction forces acting on X(t). 



velocity of the classical degrees of freedom. However, it 
was shown by Berry and Robbing that such a friction 
force is absent when the quantum system has a discrete 
spectrum. 

Several recent developments in nanoelectromechanical 
systems^-!] and spintronicJ^HUl suggest to extend these 
considerations on adiabatic reaction forces to classical 
degrees of freedom coupled to open quantum systems 
out of equilibrium. In this paper, we consider a rather 
generic model of a quantum mechanical scattering sys- 
tem (such as a coherent mesoscopic conductor within the 
Landauer-Biittiker approach EZHED see pjg [ij) which cou- 
ples to the slow classical system through the scattering 
potential. Non-equilibrium states of the quantum system 
can then be modeled by considering different distribution 
functions for the various incoming scattering channels.^ 
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Unlike the setting of Berry and Robbins, our 
scattering-theory setting naturally allows for a fric- 
tion contribution to the adiabatic reaction forces, even 
though the fast system is quantum mechanical. By the 
fluctuation-dissipation theorem, the presence of friction 
forces requires one to also include a stochastic force, and 
the classical degrees of freedom X = {X 1 , X 2 , Xjy} 
(taken to be mechanical for dcfinitcness) obey a Langevin 
dynamics, 



p cl — p 

r v — r v 



E 



>x v 



(1) 



On the left hand side (LHS), P v denotes the canonical 
momentum of coordinate X v , and we have included the 
possibility of an external classical force F cl (X). The adia- 
batic reaction forces due to the quantum environment are 
collected on the right hand side (RHS) of Eq. ([!]), where 
F(X) is the Born-Oppenheimer force exerted by the en- 
vironment, while £ denotes the stochastic Langevin force 
which represents fluctuations on top of F. The dissipative 
and Lorentz-like forces are encoded in the symmetric and 
antisymmetric parts, respectively, of the tensor 'y(X). 

The scattering approach suggests that all adiabatic 
reaction forces can be expressed in terms of the S ma- 
trix (including non-adiabatic corrections) of the quan- 
tum system. These expressions were obtained in previous 
work,™m based on a Keldysh Green's function approach 
for a closely related model. Here, our primary aim is 
to derive these expressions [given in Eqs. (48 1, |57|, (60 1, 



and (62)] directly within the setting and with methods of 



scattering theory. This alternative derivation has several 
advantages: (i) In avoiding Keldysh Green's functions 
extraneous to scattering theory, the derivation is both 
more natural and more direct, (ii) The generic scattering 
theory formulation emphasizes the generality and broad 
applicability of the results, (iii) The approach also brings 
out similarities with and differenc es fr om the seminal con- 
siderations of Berry and Robbing^! for closed quantum 
systems. 

While our model and our results are quite generic, a 
key motivation was provided by nanoelectromechanical 
system^22lMl anc i spintronic#SH2ZI devices. In these sys- 
tems, the motion of the mechanical mode or the local- 
ized spin can frequently be thought of as a slow classi- 
cal degree of freedom while the electronic conduction is 
quantum coherent and can be described as a quantum- 
mechanical scattering system, following Landauer and 
Biittiker. An important focus of recent work on adiabatic 
reaction forces in nanoelectromechanical and spintronics 
systems, often termed current-induced forces in this con- 
text, are the qualitatively new features introduced by 
out-of-equi librium q uantum environments. It is now well 
understoo d 8 * 10 * 11 *^^ that for non-equilibrium environ- 
ments, (i) the Born-Oppenheimer force is in general no 
longer conservative and thus cannot be obtained from a 
potential surface; (ii) it is possible to have negative dis- 
sipation; and (iii) a Lorentz-like force can emerge even 



for time-reversal invariant conductors. Of course, the 
approach taken here reproduces all of these results. 

Our approach may have other interesting applications. 
Since it can be applied similarly to both fcrmionic and 
bosonic environments it also provides a scattering de- 
scription of adiabatic reaction forces in optomechanical 
or cold-atom systems.^ Moreover it is also interesting to 
compare them to older results on the motion of vortices in 
superfluids. There it was shown that the geometric Berry 
phase is responsible for the Magnus force on a vortexPO 
Later it was also realized that dissipation can be obtained 
in an analogous manner within an adiabatic expansion, 
by allo wing for broadening of the energy levels of the 
systempOSS This broadening stems from the connection 
of the systems to an environment which is naturally im- 
plemented in our scattering approach. 

This manuscript is organized as follows. In Sec. [IT] 
we present the basic tools of scattering theory that are 
needed for the derivation of the adiabatic reaction forces, 
and we express the adiabatic expansion of the S-matrix 
in terms of frozen scattering states. In Sec. |III| we derive 
expressions for the adiabatic reaction forces appearing on 
the RHS of the Langevin equation in Eq. ([T]), in terms 
of the S-matrix of the quantum mechanical scattering 
system, including the first non-adiabatic correction. We 
conclude in Sec. |IV| We relegated some details to App.[A} 
and connect our results to those found in Refs.HUlandfTTI 
in App. M 



II. SCATTERING THEORY AND ADIABATIC 
EXPANSION 

A. Elements of Scattering Theory 

In this section we introduce necessary aspects of scat- 
tering theory. We consider a system described by a 
single-particle Hamiltonian H = Ho + V, where Hq is 
a free Hamiltonian and V is a scattering potential which 
depends parametrically on time through the slowly vary- 
ing classical degrees of freedom X(t). In order to de- 
scribe the system in terms of scattering states, V is as- 
sumed to be confined to a finite region in space. The 
time dependent retarded (+) and advanced (— ) scatter- 
ing states |^^(£,t)) are solutions of the time dependent 
Schrodinger equation (note that throughout this work we 
set H = 1) 



id t \v±( s ,t)) = H\y± ( e ,i)>, 



(2) 



where the index to is a combined index labeling channels 
and leads. It is convenient to define scattering states 
\ip^(e,t)) without the dynamical phase, 



IV'm(M)) 



which fulfill 



id t \ip±(E,t)) = (H-e) \i>i(s,t)) 



(3) 



(4) 
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The advanced and retarded scattering states are speci- 
fied through their boundary conditions. While the re- 
tarded state nas incoming waves only in channel to, 
the advanced state \ipm) has outgoing waves only in this 
channel, 



|^m( e >* -> Too)) = |0 m (e)) 



(5) 



where |0 m (e)) is the eigenstate of the free Hamiltonian 
in channel m, 



H \<f> m (e)) = e\4> m (e)) . 



(6) 



Eq. ([5]) holds in the weak sense, i.e., wave packets con- 
structed from the scattering states |^/^(£,t)) behave as 
free wave packets for times t — > ±00, and have energy 
e. We normalize the scattering states such that \4> m ) has 
unit flux, which implies the orthonormality relations 



<V£ (e',t)> = (0m (O) 

= 2nS (e - e') <5 mm - . 



(7) 



In the strictly adiabatic limit, the wave function 
|^(e, £)) is time-independent and hence Eq. Q reduces 
to the time-independent Schrodinger equation 



H t \i% ± {e))=e\^ ± (e)) 



(8) 



for a frozen configuration of the potential Vt = V(Kt)> 
where X t = X(t) and we have also defined H t — H(X. t ). 
We denote the frozen scattering states by |V'm* ± ( e ))- The 
superscript X t or subscript t emphasizes the parametric 
dependence on time of each quantity due to the slow 
variation of the scattering potential. Introducing frozen 
Green's functions 



Gf' A { 
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(e-H t ± irf) 



(9) 



(77 — > + ) we can write the Lippmann-Schwinger equation 
for the frozen scattering states, 

l^ie)) = \<f> m (e)) + G? /A {e) V t \<j> m (e)) , (10) 

where 10(e)) are the free eigenstates introduced in 
Eq. (|6|. Eq. (fTo]) will be of use in the next subsection. 



The frozen S-matrix St(s) is defined by the overlap 
of the frozen retarded and advanced scattering states, 
and hence depends only on the energy e of the incoming 
states, 



Sl lk {e)2TrS(e- 



Xi- 



V)K Xt+ (e)) 



(11) 



where we have isolated the singular dependence on 
energy, 8{e — e'). The frozen S-matrix is unitary, 
St(e)St(eY = 1, since scattering states are assumed to be 
normalized to unit flux. For a slowly changing system, 
the frozen S-matrix is the zeroth order contribution to 
the full S-matrix in an adiabatic expansion, as we show 
in the next subsection, and it depends parametrically on 



time through the slowly varying parameters X t . 

The exact scattering matrix for the time-dependent 
problem is defined by the overlap of the time-dependent 
scattering states introduced in Eq. ([2]), 



S nk (e',s) = (^-(e',t )\^+(e,to)). 



(12) 



The exact scattering matrix is also unitary due to the 
unit flux normalization condition 



J ^S mn {e , ,e)S nk \{e,e") = 2-K8{e'-e")8 n 



(13) 



It is important to note that the time t Q at which the over- 



lap of the scattering states is evaluated in Eq. (12), can 



be chosen arbitrarily. The independence of to can be seen 



by taking the derivative with respect to time of Eq. ( 12 ), 
and using the Schrodinger equation ([2]). This allows us to 
choose to in a convenient manner in the following section. 



We will see in Sec. Ill that even for a slow evolution, 



corrections to the adiabatic solution are important to de- 
scribe the environment-induced forces. Hence we devote 
the next subsection to calculating the first non-adiabatic 
correction to the frozen S-matrix. 



B. Adiabatic expansion and A-matrix 

The adiabatic expansion relies on the assumption that 
the classical degrees of freedom X( vary slowly in time. 
We characterize this slow time dependence by a typi- 
cal frequency f2. In finite quantum systems, adiabaticity 
requires fl to be small compared to the level spacing A. 
This condition is obviously violated in the open quantum 
systems of interest here, which have a continuous spec- 
trum. For these systems, adiabaticity requires to be 
small compared to the inverse dwell time of the electrons 
in the scattering region,^ fi <§; 1/775 . 

The adiabatic expansion is conveniently carried out in 
the Wigner representation 



S(e,t) = 



d£ 
27t' 



S(£ + £/2,e-e/2) 



(14) 



of the full S-matrix S(e',e). In the adiabatic limit, the 
S-matrix depends only slowly on the central time t. In 
fact, in the limit of a static Hamiltonian, S(e, t) becomes 
independent of t and reduces to the frozen S-matrix St (e) ■ 
For a slowly time-dependent scattering potential the 
exact S-matrix S can be writt en, up to first order in the 
adiabatic expansion, a j 10 * 40 * 41 ^ 

S(s,t) = St(e)+A t (e) + ... , (15) 

where all quantities on the RHS depend parametrically 
on time. Eq. (15) defines the A-matrix 

N 

A t (e) = '£A»(e)X v (16) 
i/=i 
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as the first-order non-adiabatic correction of S(e,t), 
which depends linearly on the velocity X and paramet- 
rically on time through X(t) . Below, we derive an im- 
portant relation which expresses the A-matrix in terms 
of the frozen scattering states \ifj' X ~ (s)}, 



Ar\e) = \{d s ^-\d v V t \^ + ) 

-\{^-\d v v t \d e ^ + ) , 



where d v = d/dX v and \d e ipf t+ ) = d e \tpf t+ ). In pre- 
vious works the A- matrix was given in terms of Green's 
function expressions j 10 l 1:L l 42 l 43 l or obtained by expa nding 
the exact solution of the time dependent problem! 40 * 41 ! 
Equation (17 1 provides a systematic way of obtaining A 



from the solution of the static scattering problem. 



To derive Eq. (17), we first compute the scatter- 
ing states \ip ± (e,t)) to first order in the adiabatic 



expansion, 



^ ± (£,i)) = |V Xt± (e)) + |^ Xt± (e)) + 



(18) 



Here, the frozen scattering state |^ Xt± ) is the zeroth- 
order term in the adiabatic expansion and corresponds 
to the strictly adiabatic limit^while |5?/' Xt± ) denotes the 
first non-adiabatic correction^. (Here, we omit the chan- 
nel index for notational simplicity.) Inserting Eq. (18) 
into the Schrodinger equation Eq. Q , using Eq. Q , and 
comparing terms of first order in the adiabatic expansion, 
we find 



id t \^ ± {e)) = {H t -e)\5^ ± {e)) 



(19) 



where dt indicates the parametric derivative with respect 
to time. With Eq. we conclude that \6->p Xt± ) = 



-iG?^dt\^ 

Eq. @, 



X; ^ and hence, plugging this back into 

\^{e, t)) = |^ Xt± (£)) - i G? /A {e) d t |V Xt± (e)) + • • • • 

(20) 

We can also express 9t|-0 Xt± ) in terms of the frozen 
scattering states by using the Lippmann-Schwinger equa- 
tion given in Eq. (10). Taking the time derivative of 



r,R/A 



R/As /r ,R/A 



= Gt'VtG 



where 



Eq. (10) and using that G t 

Vt — dtVt, the time derivative of the scattering states 



can be expressed as^ 

9 t |V Xt± (e))=Gf M ( £ )\> t |V Xt± (£)) 



(21) 



and hence we obtain the desired result from Eq. (20), 

\^(e,t)) = 1^(6)) - i (Gf^V^ie)) + ..., 

(22) 

which is valid to first order in the adiabatic expansion. 

With this expansion of the scattering states, the adia- 
batic expansion of the full S-matrix can be performed, 



starting from the definition Eq. ( 12 ) and the Wigner 



transform given in Eq. ( 14 ) . It is convenient to use the 
arbitrariness of t in Eq. ( 12 ) by choosing i as the cen- 



tral time, to = t. Defining e± = e ± e/2, the Wigner 
transformed S-matrix can now be approximated to first 
order in the adiabatic expansion as 



(17) S(e,t) = 



(*- (e+,t)|* + (£-,*)) 
^ Xt -(e + )|^ x '+(e_)) 



de 
— < 
2tt 

27r' 

g(^-( £+ )|[Gf( £ _)] 2 T4|V^ + ( £ -)) 
-,-//g(^-( £+ )|]/ t [Gf( £+ )] 2 |^ + (e-)) 



(23) 



where we have used Eq. ( 22 ) . We now employ the iden- 
tities Gf (e ± )\tp Xt + (e T )) = {±e + irf}- 1 ^^ (e T )) and 
[(£+ ir))~ 2 - (— e + iri)~ 2 ] = 2Tridg5(e) to obtain 

S(e,t) = S t (e) 

de [8,6(6)} (e+) \V t \^ + (e_) ) + . . . . 

(24) 



Integrating by parts with respect to e, we find 
S(e,t) = S t (e) + \(d e ^- (e) \V t \^+ (e) ) 
- l(^'(e)\Vt\d e ^ + (e)> + ... , 



(25) 



which gives the full S-matrix S in terms of the frozen 
S-matrix St defined in Eq. (11), and the first non- 
adiabatic correction matrix A t (A-matrix) as anticipated 
in Eq. (fl7}P 



We finish this section by deriving some identities for 
the S and A-matrices that will be of use in the derivation 
of the adiabatic reaction forces. The frozen S-matrix can 
be written as 

Sf{e) = S nk ~ i(4fr-{e)\ V t \Me)) , (26) 

which follows from recognizing that the second term on 
the RHS ofj26} is the frozen T-matrixP This together 
with Eqs. (10) and (21) gives the time derivative of the 



frozen S-matrix in terms of the frozen scattering states, 
d t Sf(e) = -t(^-(e)\V t \^f t+ (e)) . (27) 

Hence 

8 e (^- (e) |T4|V Xt+ (e) ) = id £ d t S t (e) (28) 
and we obtain an alternative expression for the A-matrix 
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by comparing with Eq. (17 1 



A t {e) = -<^ x - (e) \V t \d £ ^ + (e) ) + -d e d t S t {e) . 

(29) 

The S and A-ma trices are related through unitarity of 
the exact S-matrixj 10 ! 11 * 40 ^ resulting in the identity 



S\A t + A\S t = l - (d t S\d e S t - d e Std t S t 



(30) 



where all quantities are evaluated at the same energy. We 
can check that our explicit expression for the A-matrix 



in Eq. (17 1 indeed fulfills this condition. This is shown 
in App.|£J 



corresponding basis. Hence we have 

Ht = /£/ %T.^nk^t)a k {e',t). (33) 

mk 

It is straightforward to show that the quantity 
(aj n (e, t)a k (e' , i)) is independent of time, by noting that 
both the retarded scattering states \^^(e,t)) and the 
density matrix p evolved unitarily from the unperturbed 
states. Then, the occupation f n (e) of a scattering state 
in channel n is governed by the corresponding reservoir, 
as in the Landauer-Buttiker theory of mesoscopic con- 
ductors, 



(alnie, t)a k (e', t)) = f k (e)5 km 2ir5(e - e') . 



(34) 



III. ADIABATIC REACTION FORCES 

A. Adiabatic reaction forces and scattering states 

The force operator in the Schrodinger picture can be 
defined as 



•Fx = -V7Z . 



(31) 



where the gradient is taken with respect to X and H = 
H(X) is the (non-interacting) many-body Hamiltonian 
of the quantum system. Notice that H includes terms 
arising from the free Hamiltonian Hq of the fast degrees 
of freedom and the scattering potential V(X) which de- 
pends parametrically on the slow, classical variables X. 

Then for a given trajectory X t the average force that 
the out-of-equilibrium quantum environment exerts on 
the classical degrees of freedom X at time t is given by 



jr(t) = j-[x t ] = <Jx t ; 



(32) 



Here the expectation value indicates quantum-statistical 
averaging for a given trajectory X t : (...) = Tr{p(t) . . .}, 
where p(t) is the many-body density matrix of the sys- 
tem at time t. Notice that p{t) and thus the force T are 
junctionals of the trajectory X t , and therefore J 7 depends 
on time through X t and its time derivatives, X t , X t , .... 
In fact, the adiabatic expansion consists of making a sys- 
tematic expansion in these latter quantities. Also, equa- 
tion (32) gives only the average force. As mentioned 



in the Introduction, the Langevin dynamics includes the 
stochastic fluctuations of — VH, which we will consider 
further below. 

To compute the quantum-statistical average (. . .) in 
Eq. ( 32 ) , we write the many-body Hamiltonian in 
terms of creation and annihilation operators a^(£,i) and 
a n (e,t) that create/annihilate the retarded scattering 
states \^>^(e,t)). Since the time evolution is unitary, the 
retarded scattering states constitute a complete basis at 
any time t. Note that we are working in the Schrodinger 
representation, and the time t appears as a label in the 
creation/annihilation operators a simply to identify the 



Expressing the force operator J-"x t = —WHt in terms of 
these creation and annihilation operators as 

VHt = /£/ ^T,FH t ) mk al(s,t)a k ( S >,t), (35) 



we are now in a position to evaluate the average adiabatic 
reaction force as 



T 



E 

k 



d£ 

2^ 



/ fe ( £ )(^+(M)|Vff t |V<+(M)). (36) 



This expression allows us to perform an adiabatic expan- 
sion of the reaction force using the adiabatic expansion 
of the scattering states developed in Sec. |TTj 

Inserting the adiabatic expansion of the scattering 
states given by Eq. (22) into Eq. (36), and keeping terms 



up to first order in the adiabatic expansion, we find 
de 



JX f = - 



2tt 



f k (e)(^ + (e)\VV t \^ + (s)) 



-i J ^fk(s){^ + (e)\d t V t (Gt) 2 VV t \^f' + (e)) 
J ^Ue)(^ + {e)\VV t (G?) 2 d t V t \^+{s)) 



+ i 



(37) 

where we have used that VH t = VVt and left the sum- 
mation over the channel index k implicit. 



Equation (37) yields the deterministic reaction forces 
appearing on the RHS of the Langevin equation, Eq. (|l 
The zeroth order term, given by the first line of Eq. (137 
is independent of the velocity X and gives the (possibly 
non-conservative) Born-Oppenheimer force F(X). The 
first-order contribution, given by the second and third 
terms of Eq. (37), represents the forces that depend lin- 

7 (X)-X. 
To define 



early on the velocity of the classical modes, - 

Let us now turn to the force fluctuations 
time-dependent force fluctuations at the quantum me- 
chanical level, one needs to go to the Heisenberg picture, 
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J-x — > ^x(i) and define the Heisenberg force fluctuation 
operators, 



(38) 



There are two different contributions to the stochastic 
force: (i) fluctuations at finite temperatures and (ii) non- 
equilibrium noise which is a consequence of the proba- 
bilistic nature of the scattering process. Since the quan- 
tum degrees of freedom are fast compared to the mechan- 
ical motion, the correlator D a p(t,t') of the stochastic 
force is local on the relevant time scales of the Langevin 
equation ([I]), 

D a0 (t,t') = {(Ut)£ p (t'))} s ~D aP (t)S(t-t f ), (39) 

were the subscript s denotes symmetrizing with respect 
to a and (3. To account for these fluctuations and to sat- 
isfy the fluctuation-dissipation theorem for the classical 
variable X, one must incorporate in Eq. ([!]) the cla s- 
sical stochastic force terms, £ a {t), obeying ^ a (t)^p(t') = 
D a p(t)5(t—t') where the overline corresponds to the clas- 
sical averaging implicit in the Langevin equation. . 

In order to determine the correlator D a p (t) , we average 
D a p(t,t') over the fast degrees of freedom corresponding 
to the relative time r, 



D aP (X t ) = / drD 



(40) 



It is sufficient to evaluate this correlator in the fully adia- 
batic limit since this already ens ures t hat the fluctuation- 
dissipation theorem be satisfied! 10 ^ 11 ! Hence we freeze the 
value of X t (and by that the Hamiltonian) , and evaluate 
the force fluctuations with this static Hamiltonian. Then 
we can work in the frozen scattering state basis |t/> x ), 
where the Schrodinger force operator can be expressed as 

^ = ~ / £ / t(£) [daHt] - k ■ (41) 

Here the channel summation is again implicit and a^*^ (e) 
creates a frozen (retarded) scattering state in channel 
m with energy e at time t, i.e., an eigenstate of the 
frozen Hamiltonian H t . The superscript X t indicates 
the strictly adiabatic condition. Within this adiabatic 
approximation, changing to the Heisenberg picture sim- 



ply amounts to replacing in Eq. (41) the operators a 



by the corresponding Heisenberg operators, a k ' (s) — ¥ 

a*<(M + T) = c- ier a^(£)PThe correlator D a0 (t) can 
then be calculated by means of the identity^ 



(al n (e 1 )a n (e 2 )al{s 3 )a l {e4))-(al l {e 1 )a n (e 2 ))(al(e 3 )ai(e 4 )) 

= (2tt) 2 / m (ei) [I T fk{£2)] 5 m i S nk S(ei - e 4 ) S(e 2 - e 3 ) , 

(42) 

where the upper sign (— ) refers to fermions and the lower 
sign (+) to bosons, and for simplicity, the label X t has 



been dropped. Applying Eq. ( 42 1 to Eq. ( 40 1 we finally 
arrive at 

D af} (X t ) = JdrJ ~/ m (e) [IT /*(£)] 

x { <V^ t+ (e) | d a H t |^ Xt+ (e)) <V> Xt+ (e) | d p H t \ ^ + (e))} 

(43) f 

where indicates symmetrization with respect to 

the indices a, (3. In the remainder of this section, we shall 



express the average adiabatic reaction force in Eq. (37) 



as well as the correlator in Eq. ( 43 ) in terms of the S and 
A-matrices. 



B. Born-Oppenheimer force 

The Born-Oppcnhcimer force F can be expressed solely 
in terms of frozen scattering states. From Eq. (37), it is 
given by 



F„ = - 



de 
27 



i-i (44) 
For turning Eq. (44) into an expression involving the 



S-matrix, we insert a resolution of the identity 1 = 

Fa = - ft I %Y.^){^He)\^-(e')) 

x(^-(e')\d a V t \^+(e)). 

(45) 

Since the dependence on time t is parametric through 
X t , from Eq. ((27|) we have 



(46) 



d a S nk (e) = -i(^~(s)\d a V\^f t+ (e)) 



Putting this together with the expression for the frozen 
S-matrix Eq. (11) we obtain 



2ni 



J2fn(e)S} nm (e)d a Sr(e) 



(47) 



In matrix notation, this gives the Born-Oppenheimer 
force 

F a (X t )= [ ^.J2fn(e)tr{n n Sl(e)d a S t (e)} , (48) 



where tr{. . .} denotes a trace over scattering channels, 
and n„ is a projector onto channel n. Eq. (48) coincides 



with the one obtained in Ref. \W\ via a non-equilibrium 
Keldysh calculation for the current-induced-forces in a 
nanoelectromechanical system. 



The expression given in Eq. ( 48 1 can be motivated 
by realizing its connection with the Friedel sum rule.^ 



Considering a finite system with discrete energy levels for the Friedel sum rule in terms of S-matrices; 
El — i? 4 (X t ), the Born-Oppenheimer force in equilib- 
rium is given by 
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F a (X t ) = -J2f{E l t )d a El 



(49) 



If we replace E\ -> / deeS(e - E\) in Eq. ([49]), 



we can take the limit of the system size to infinity 
by writing the number of states up to energy e as 
N(e, Xj) = /f deV(e', Xt) where f is the density of 
states and we have used the identity d a Q (e — El) = 
-S (e - Ej) d a Ej. In this limit, Eq. d49l) takes the form 



F a (Xt)= / def(e)d a N(e,X t ). 



(50) 



The quantity d Q A~ is known as the emissivitjBS and 
plays a key role in the problem of adiabatic quantum 
pumping Psl [Note that expressions of the type Sld a S t 
also appear in the context of quantum pumping as "re- 
sponse matrices" .] Making use of the general expression 



N(e,X t ) = -L t r{lnS 4 ( £ )} 
2m 



(51) 



the emissivity can be expressed as 

d a N(e,X t ) = ^-tr{sl(e)d a S t (e)} , (52) 



and hence we recover Eq. (48) for the particular case 



that the system is in equilibrium. 



C. Friction and geometric magnetic force 

We now express the tensor 7 in terms of the frozen 
S-matrix and the first order non-adiabatic correction, 
the A-matrix. The first order correction to the Born- 
Oppenheimer force is given by the two last lines of 
Eq. (37). With the aid of the chain rule d t = X a d a , 



it is straightforward to show that 



l aP = ij ^f k (e)(^ + \d fi V t (Gt) 2 d a V t \^+)-i J ^f k (s)(i>f t+ \d a V t (G?) 2 dpV t \i>ft + ), (53) 

where we have omitted energy variables and left the sum over k implicit. We split the tensor 7 into a symmetric 
part j^b — 1/2 (7^ + 7s a ), corresponding to the friction force, and an antisymmetric part j® g — 1/2 (jap — Ipa), 
corresponding to the emergent Lorentz force. 

We first consider the symmetric, dissipative contribution 

Y a p=i J ^fk(e)(^ + \d a V t (G?) 2 d p V t \^ + ) s -z J ^-f k{ e)(^ + \d a V t (G?) 2 dpV t \^ + ) s , (54) 
Using the identity (Gj 4 ) 2 — (Gf-) 2 = —d £ (Gf — Gj* ) = —2i:id e 5(£ — H t ) and integrating by parts we obtain 

Y a p= [ &e[-d £ f k {e)](^ + \d a V t 5(E ~ H t )dpV t \^ + ) s - [ def k (e)(d^f t+ \d a V t S(e-H t )d^V t \^ + ) s 

J (55) 

def k (e)(^f t+ \d a V t S(e - H t )d p V t \d e ^ + ) s . 

We now insert a resolution of the identity J^i J i^lV'i Xt ( £ '))(V'; Xt_ (£')l between the two potential terms and find 

de 



lU = ~ J i d ^ £ ) {<V )Xt+ |9^t|^ Xt ->(^ Xt "|^F t |^ + >} s 

-f k {e) Ud e ^ + \d a V t \^~)(^-\d V t \^ + )\ (56) 

7T ^ ) s 



ds 

27 

de 



^fk(e){(^ + \djr t \^-)(^-\d fi V t \a e ^*+)} 



Comparing this expression with Eqs. (|29[) and (|27[) and using the definition of the A-matrix given in Eq. (16), leads 



immediately to the final result 

7^(X t ) = J g^[-9 e /„( e )]tr{n„9 a 5l( £ )^5 t ( £ )} s 



(57) 



Equation ( [57] ) recovers the frictional force obtained first 
in Ref. [TOJ Thus, we conclude that the classical de- 
grees of freedom are indeed subject to a friction force 
due to the coupling to a quantum mechanical scattering 
system. This is in stark contrast with the coupling to 
a finite quantum system where Berry and Robbins find 
that the frictional contribution to the adiabatic reaction 
force vanishes!^ 

The antisymmetric part of the damping matrix has 
the role of an effective orbital magnetic field acting on 
the multidimensional space of X. From Eq. (53), it is 
given by 



7a/3 



f(e)(^ + \d a V t (G?) 2 d p V t \^ + ) a 



de 
2ttI 



(58) 

To evaluate this expression we observe that, using 
Eq. ((29]) and (Gf ) 2 = -d e Gf, it follows straightfor- 
wardly that 

[d a A f3 t ) a = {lF*-\d p V t {G?fd a V t \iP x < + ) a , (59) 
with an analogous expression involving Gf. Similar ma- 



nipulations to the ones employed in Eq. ( 56 ) lead to the 
result for the Lorentz-like term of the -y matrix 

7 ^ (Xt) = / £ £ Ue) 

J n 

• {n„ [sfe)df,A*(e) - d p Af(e)S t (e)] . 

(60) 



x tr • 



This expression agrees with the one obtained in Ref. 1101 
Note that within the scattering formalism all the above 
relations are well defined and we do not encounter any 
divergences in contrast to Ref. [HI The reason for this 
is that the particles spend a finite time in the scattering 
region, which is implicit in the time-dependent scattering 
formalism. 



D. Stochastic force 

The stochastic force can be written in terms of the 
frozen S-matrix by inserting two resolutions of the iden- 
tity of the form 1 = J Jf £ fc ^"(e))^ into 
Eq. (43). We can then identify the frozen S-matrix by 



use of Eq. (11) as well as its derivatives given by Eq. ( 26 1 . 
This yields 

D a0 = [ p-f n (e) [lT/™(s)] 

7 "" (61) 
X ^d a St nk {e)S k r{e)St ml {e)d p S\ n {e) 

In matrix notation, this can equivalently be written as 

Atf(Xt) = £ J g/„( £ ) [it f m (s)} 

nm 

x tr |n„ [S\{e) d a St(e)] ^ U m Sj(e) d p S t {e)\ , 

(62) 

which agrees with the expression in Ref. UOI when dealing 
with fermions. It can be shown that this expression ful- 
fills the fluctuation-dissipation theorem in equilibrium,^ 



D = 2T~f s , with 7 S given by Eq. (57) and evaluated in 
equilibrium. 



IV. CONCLUSION 

Slow degrees of freedom coupled to a fast quantum sys- 
tem are subject to adiabatic reaction forces. Currently, 
these forces play a pivotal role in the context of nano- 
electromechanical and spintronics systems with a slow 
mechanical mode or spin coupled to fast electronic de- 
grees of freedom. While early work on the adiabatic re- 
action forces focused on closed quantum systems, nano- 
electromechanicals and spintronics typically involve open 
electronic systems driven out of equilibrium by voltage 
sources. 

These developments have motivated us to consider adi- 
abatic reaction forces for a generic model of a slow clas- 
sical degree of freedom coupled to a quantum mechan- 
ical scattering system. Non-equilibrium is incorporated 
into this (non-interacting) many-body model by assum- 
ing that the filling of the incoming scattering channels 
is controlled by various reservoirs. In the context of na- 
noelectromcchanics and spintronics, this model follows 
naturally if the electronic degrees of freedom take the 
form of a mesoscopic Landauer-Biittiker conductor. 

It was sho wn rec ently within a Keldysh Green's func- 
tion approac h 1 10 * 11 ! that the adiabatic reaction forces can 
be expressed entirely in terms of the adiabatic S-matrix 
and its first non-adiabatic correction, the A-matrix. The 
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main result of the present paper is an alternative deriva- 
tion of these results within the setting and with the meth- 
ods of scattering theory. 

In addition to being more natural and more direct, this 
derivation has several further advantages. To start with, 
we present useful expressions for the A-matrix in terms 
of the adiabatic scattering states which should simplify 
its calculation for specific applications. The general set- 
ting within the context of scattering theory facilitates 
comparison with the earlier results on adiabatic reaction 
forces for closed quantum systems. Most prominently, 
there is no frictional force for closed quantum systems 
while such a force emerges naturally for a quantum- 
mechanical scattering system. Moreover, the approach 
clarifies the limits of validity. While for closed quantum 
systems, the adiabatic condition involves the level spac- 
ing, the latter is replaced here by the dwell time of the 
fast system in the scattering region. Finally, the general 
setting emphasizes the generality and wide applicability 
of our results. The fast quantum system can be fermionic 
as in nanoclcctromcchanics, a spin degree of freedom as 
in spintronics, or bosonic as in optomechanics or cold- 
atom systems. 
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Appendix A RELATION BETWEEN THE S AND 
A-MATRIX 



Here we prove the identity given in Eq. ( 30 ) . For better 



readability we omit the channel index and energy depen- 
dence (all quantities are evaluated at the same energy e). 
Starting with Eq. ( 29 1 we find 



St A t + A\S t 



-d e {^ + \v t \^ + ) + l - [s\d e d t s t - d £ d t s\s t 



S\dtSt — dtSjSt 



d t sld e s t - d e s\d t s t 



s\d e d t s t - d e d t s\s t 



ity in Eq. (A.l). 



where we made use of Eq. ( 27 ) to obtain the second equal- 



Appendix B APPLICATION: A QUANTUM DOT 
COUPLED TO LEADS 



As we mentioned above, the S and A-matrix expres- 
sions presented in this work for the reaction forces were 
obtained first, within a different formalism, in Ref. 1 101 for 
the forces that an applied current exerts over the slow 
vibrational degrees of freedom of a nanomechanical oscil- 
lator. In this section we show how our formalism relates 
to the one presented in Ref. QIT| where the Hamiltonian 



H(X) = H x +H L + H L 



(B.l) 



was considered, which models a quantum dot connected 
to leads. The "heavy" classical degrees of freedom X(i) = 
{Xx{t),X%(t) . . .Xjv(t)} in this case are the mechanical 
vibrational modes of the dot, which couple to the elec- 
trons in the dot. The different terms of the Hamiltonian 



in Eq. (B.l I are given by 



p2 

2M„ 



U(X) 



H D = J24n IMX)] W d m , 

mm' 

H T = J -~^^(c\ ] {e)W vm {e)d m + h.c) 



(B.2) 

(B.3) 
(B.4) 

(B.5) 



Hl models the leads, where cL(e) [cr,(e)\ creates [anni- 
hilates] electrons in a flux normalized state \4> n (s)) in- 
coming from rj with chemical potential ji^ (77 combines 
channel and lead index, the chemical potential depends 
only on the lead index). Hx represents the free evo- 
lution of the mechanical degrees of freedom of the dot. 
He, is the Hamiltonian of the dot, containing the elec- 
tronic levels plus the coupling of the electrons in the 
dot to X via a general function /i(X). The operators 
d m (d m ) create (annihilate) a dot-electron in the state 
|m). Finally, Hx indicates the tunneling process between 
the dot's levels and the leads with tunneling amplitude 



W vm (e) = (<£„(£) I W|m>/V27r. 



The electronic part of the Hamiltonian in Eq. (B.l) 



can be interpreted as a scattering problem, where the free 
Hamiltonian is given by Hl and the dot defines a scat- 
tering potential V = U D W^U L + Il L Wn D + Il D H D U D 
where n^, and project onto the lead and dot space, 
respectively (Note that • TId = IId • Hi, = 0.) 

We can then write the Lippmann-Schwinger equation 



as 



(A.l) |< t + (e)) =n L |^(e))+Gf(e)y t n L |0 r ,(e)> (B.6) 



where Gf'(e) — (e — H t + irf) 1 is the frozen Green's 
function of the dot plus lead. The projection of Eq. (B.6 ) 
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onto the dot space takes the form 

U D K t+ (e)) = n D Gf (e) U D Wl U L \<j> v (e)) 

= G^(s)W^\^(e)) (B.7) 

with Gp(e) =Ud Gf (e) Hd, the dot's frozen Green func- 
tion. 

We are interested in an explicit expression for the 
A-matrix. Since d t V t — ^DdfHjjIljj, we obtain from 



Eq. (17) 



Using Eq. ( |B.7[ ) this can be cast into the form 

AT = n{d E (W vk [Ggy [d t h lm ] [G R D ] mn Wl 
-W vk [G B D ] kl [dth lm ]d s ([G%] mn Wl)} , 



(B.9) 



where the summation over repeated indices is implied. 
Applying our formalism to the Hamiltonian (B.l) there- 



fore indeed recovers the expression of the A-matrix given 
by Ref . |TJB 



ATie) = 2 [(d £ ^-(e)\ll D d t H D U D \^+(s)) 
-«* _ (e)| Tl D d t H D U D \d e ^ + {e))\ 



(B. 
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